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Abstract. We use far-UV absorption spectra obtained with FUSE towards three late B stars to study the formation 
and excitation of H2 in the diffuse ISM. The data interpretation relies on a model of the chemical and thermal 
balance in photon-illuminated gas. The data constrain well the n R product between gas density and H2 formation 
rate on dust grains: n R = 1 to 2.210~ 15 s~\ For each line of sight the mean effective H2 density n, assumed 
uniform, is obtained by the best fit of the model to the observed N(J=1)/N(J=0) ratio, since the radiation field 
is known. Combining n with the nR values, we find similar H2 formation rates for the three stars of about 
R = 410 _17 cm 3 s _1 . 

Because the target stars do not interact with the absorbing matter we can show that the H2 excitation in the 
J > 2 levels cannot be accounted for by the UV pumping of the cold H2 but implies collisional excitation in 
regions where the gas is much warmer. The existence of warm H2 is corroborated by the fact that the star with 
the largest column density of CH + has the largest amount of warm H2. 
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1. Introduction 

The H2 formation is a key process for the understanding of 
the thermal and density structure as well as the chemical 
evolution of the interstellar medium (ISM) . The H2 forma- 
tion rate was first estimated through the modelling of the 
hydrogen recombination on dust surfaces (e.g. Hollenbach 
et al. 1971). Based on an analysis of Copernicus ob- 
servations of atomic and molecular hydrogen in the lo- 
cal diffuse clouds, Jura (1975a) proposed an H 2 forma- 
tion rate (R = 3 10 _17 cm 3 s _1 ) which corresponds to the 
Hollenbach et al. (1971) prediction for a total grain surface 
of 10 _21 cm 2 /H and a recombination efficiency of 0.5. 

The excitation of the H2 rotational levels from the 
ground state observed in absorption in the UV is a di- 
agnostic of physical conditions. In diffuse clouds, the low 
J lines provide a measure of the gas temperature while the 
excitation of the J > 2 levels is generally interpreted as 
a result of the fluorescence cascade following H2 pumping 
by the UV radiation from the target stars. But collisional 
excitation in shocks driven by the star have also been con- 
sidered. 

Since the pioneering work of Black & Dalgarno (1976) 
on which the Jura analysis of Copernicus data is based, 
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much progress has been made in the modelling of H2 in 
space, in particular about the fluorescence cascade after 
UV pumping of electronic transitions and collisional de- 
excitation rates (Combes & Pineau des Forets 2000). The 
Far Ultraviolet Spectroscopic Explorer (FUSE) is also now 
providing new UV absorption observations of Galactic H2 
superseding the Copernicus observations by their sensi- 
tivity (Snow et al. 2000, Shull et al. 2000, Rachford et 
al. 2001). The topic of H2 formation and excitation has 
also been revived by the observation of the mid-infrared 
transitions between the rotational levels of the vibrational 
ground state. These data have been used to estimate the 
H2 formation rate in warm photo-dissociation regions at 
the surface of molecular clouds (Draine & Bertoldi 1999, 
Habart et al. 2002) where the gas and the dust are both 
warmer than in the diffuse ISM and thus where the H 2 
formation efficiency or the processes involved might dif- 
fer. ISO observations are not sensitive enough to detect 
the mid-IR H2 line emission in the low to moderate col- 
umn density lines of sight studied in the UV but they have 
allowed to detect an extended warm H2 component away 
from star forming regions across the Galaxy (Verstraete et 
al. 1999) and in the edge-on galaxy NGC 891 (Valentijn 
& van der Werf 1999). The Galactic data has been inter- 
preted as evidence for the existence of warm H 2 gas heated 
by the dissipation of kinetic turbulent energy. 



2 



C. Gry et al.: H2 formation and excitation in diffuse fSM 



In this paper, we re-consider the question of H2 for- 
mation and excitation in the diffuse ISM by analysing 
FUSE observations of three late B stars located behind 
the Chamaeleon clouds. IRAS images show that these 
stars unlike most earlier type stars usually observed in the 
UV are truely background field stars that do not interact 
with the matter responsible for the absorption and that 
do not contribute to the incident radiation field. These 
three stars were part of a larger sample of Chamaeleon 
lines of sight observed with IUE to correlate the UV ex- 
tinction curve with changes in the dust size distribution 
traced by the IRAS mid-IR to far-IR colors (Boulanger, 
Prevot, Gry 1994). The interstellar medium along these 
lines of sight was further characterized with high reso- 
lution optical and ultraviolet absorption spectra with a 
spectral resolution of 10 5 obtained at ESO and UV spec- 
tra with the Goddard high resolution spectrograph on the 
Hubble Space Telescope (Gry et al. 1998). Characteristics 
of these lines of sights and previously observed column 
densities are given in Table @. 

The FUSE data results (Sect. ||) are analysed by 
modelling the interaction of the gas with the Solar 
Neighborhood radiation field. This analysis leads on one 
hand to an estimate of the H2 formation rate (Sect. ||) 
and on the other hand to a discussion of the H2 excita- 
tion (Sect. ^) where we argue that H2 column densities at 
J > 2 cannot be acounted for by UV pumping. 




1090.00 1091.00 1092.00 1093.00 1094.00 1095.00 1096.00 
wavelength 



Fig. 1. Example of fit of the H2 J=0 and J=l level tran- 
sitions in the line of sight towards HD102065. The plot 
shows the FUSE observations, the synthetic absorption 
components due to the individual transitions, and the re- 
sulting profile. Only one band (Lyman(l,0)) is shown but 
all available bands are fitted together. Only the column 
density of the J=0 and J=l levels (marked H20 and H21) 
are derived from the fit. The higher H2 levels (in this case 
J=2, J=5 and J=6) as well as the HD lines are included in 
the fit, but their column densities cannot be derived from 
this fit. 



2. Observations and data analysis 

High resolution far UV (905-1187 A) spectra of 
HD102065, HD108927 and HD96675 have been obtained 
in April and May 2000 by FUSE observations of re- 
spectively 6.8, 6.5 and 10.2 ks, as part of the PI Team 
Guaranteed Time. FUSE design, data processing and per- 
formance have been described in Moos et al. (2000) and 
Sahnow et al. (2000). The data have been reprocessed with 
the calibration pipeline version 1.8.7. 

To derive the H2 column densities, we have used the 
bands located in the best part of the spectra, namely the 
Lyman (0,0), (1,0), (2,0), (3,0), (4,0) and (5,0) bands. In 
effect the UV flux of these stars decreases rapidly with 
wavelength shortward of Ly/3, due to their relatively late 
types and their interstellar extinction. The stellar Ly/3 ab- 
sorption is so broad and deep in these spectra that it com- 
pletely wipes away the interstellar line, prohibiting the 
determination of HI column densities. 

The J=0 and J=l H2 lines are heavily saturated and 
their Lorentzian wings allowed us to derive precise column 
densities for these two levels independently of the gas ve- 
locity distribution. They have been measured through pro- 
file fitting with the program "Owens" which allows to fit 
together (blended) lines from different H2 levels and other 
species. All available bands are fitted together meaning 
that the derived column density gives the best fit on all 
bands together. This allows to decrease the errors mainly 
due to the uncertainty on the continuum placement be- 
cause of overlapping wings of adjacent bands. The error 



bars on column densities have been estimated empirically 
as the dispersion of the results of five fits performed with 
different sets of four H2 bands each. 

The derived H2 column densities are listed in Table fn. 
An example of the data and the fits is shown in Fig.l^ 
for the Lyman(l,0) band in the line of sight towards 
HD102065. 

For the higher excited levels, the lines are less satu- 
rated and fall on the flat part of the curve of growth where 
the column densities depend heavily on the gas velocity 
distribution. Thanks to high resolution optical observa- 
tions of the molecules CH and CH + , we have derived the 
b-value of the gas responsible for the CH absorption for 
the three lines of sight and of the gas responsible for the 
CH+ absorption for two of them (Table 0). For HD102065 
we also have a measurement of the C 1 b-value from new 
high resolution HST-STIS data: b(Ci) = 1.8 ±0.1 km/s, 
close to b(CH). We derive the high-J H2 column densities 
via a curve of growth analysis by adopting successively the 
CH and CH + b- values. The results are given in Table || 
and the method is illustrated in Fig. || in the case of the 
line of sight toward HD 102065 for which a significantly 
different b- values have been measured for CH and CH + . 
Fig. |^ shows the fit of all measurements for levels 3—2 to 
J=5 to the curve of growth corresponding to b(CH) (b 
= 1.9 km" 1 , top) and b(CH+) (b = 3.0 km -1 , bottom). 
One must keep in mind that significant uncertainties can 
be attached to these column densities in case the velocity 
distribution of the excited H2 gas is different from that of 
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HD102065 


HD108927 


HD96675 


badopted 




1.9 km/s 


3.0 km/s 


2.2 km/s 


1.6 km/s 2.2 km/s 


N(H 2 ,J= 


=2) 


2.6 10 18 


2.5 10 18 


1.7 10 17 


3.3 10 18 3.4 10 18 


N(H 2 ,J= 


=3) 


3.1 10 17 


1.1 10 17 


9.3 10 16 


1.7 10 17 8.7 10 16 


N(H 2 ,J= 


=4) 


5.610 1B 


6.0 10 15 


1.5 10 15 


4.0 10 15 1.4 10 15 


N(H 2 ,J= 


=5) 


1.2 10 15 


4.5 10 14 


< 6 10 13 


< 2 10 14 



Table 2. H 2 column densities (in cm~ 2 ) in the excited (J>2) levels, derived from a curve of growth analysis by 
assuming successively b = bcH and b = b CH +, when available. It is likely that the true N(H2) is intermediary or close 
to one of these two estimates. 





HD102065 


HD108927 


HD96675 


Sp. Type 


B9IV 


B5V 


B7V 


d(pc) 


170 


390 


250 


E(B-V) 


0.17 ± 0.04 


0.23 ± 0.04 


0.31 ± 0.03 


A v 


0.67 ± 0.12 


0.68 ± 0.10 


1.1 ± 0.15 


Rv 


3.9 ± 0.4 


3.0 ± 0.2 


3.5 ± 0.2 


N(CH) 


6.3 10 12 


1.4 10 13 


2.2 10 13 


b(CH) 


1.9 km s _1 


2.2 km s" 1 


1.6 km s _1 


N(CH+) 


1.2 10 13 




2.8 10 12 


b(CH+) 


3.0 km s" 1 




2.2 km s" 1 


N(CO) 


(7 ± 3) 10 13 




> 10 15 


N(Ci) 


(4 ± 1) 10 14 




> 1.7 10 15 


Ntotal 


9.9 10 20 


1.3 10 21 


1.8 10 21 


N(H 2 ,J=0) 
N(H 2 ,J=1) 

f _2N(H 2 ) 
J N(H) 


2.0±0.2 10 20 
1.4 ±0.1 10 20 

0.69±0.12 


2.0±0.2 10 20 
1.2±0.110 20 

0.49±0.09 


4.9±0.910 20 
2.1 ±0.6 10 20 

0.76 ±0.15 



Table 1. Interstellar absorption toward the stars. 
Column densities are in cm~ 2 . 

CH and CH + column densities were derived from high 
resolution (R=110 000) optical spectra obtained with the 
ESO 3.6 m telescope and the Coude Echelle Spectrometer. 
The C 1 column densities were derived from HST GHRS 
observations (Gry et al. 1998). 

The total column density, Ntotal , has been derived from 
E(B-V), N to t a i = 5.8 10 21 E(B-V). 



CH and CH + . However in view of the general good linear 
correlation of N(CH) with N(H 2 ) at these column densi- 
ties (Mattila, 1986), and the correlations of N(CH + ) with 
excited H 2 column densities (see Sect. [I]), it is likely that 
the J>2 column densities are intermediary or close to one 
of these two estimates. 

3. H2 formation 

3.1. Determination of the nR product 

The local balance between H 2 formation and photo- 
dissociation can be written as: njji nR — n^ 2 Po G S, 
where njji, i-h 2 ? n are the atomic, molecular and total 
hydrogen densities (n = rijji ± 2 njj 2 ) , R the H 2 forma- 
tion rate, /3q is the Solar Neighborhood value of the H 2 
photo-dissociation rate in the absence of shielding, G is 
the radiation field value in Solar neighborhood units and 
S a shielding factor including dust extinction and H 2 self- 



Curves of Growth for HD102065 




Fig. 2. Curves of growth for the excited levels of H 2 . The 
top plot shows the curve of growth for b(H 2 ) = b(CH) 
and the bottom plot that for b(fi£ = b(CH+). The J> 1 
column densities listed in Table |lj have been derived by 
fitting the measurements of all J> 1 levels together to 
both curves of growth respectively. 

shielding. For the stars studied here, the IRAS maps show 
that the stars do not heat the matter responsible for the 
absorption, we can thus assume that the stars do not 
contribute to the radiation field and we can take G=l 
(Boulanger et al. 1994). For a constant n, integration over 
the line of sight leads to: nR = < S > where 

/ is the molecular hydrogen fraction: f=2 N(H 2 )/N to tai, 
where N to tai, the total column density is derived from the 
extinction, N tota i = 5.8 10 21 E(B-V), and < S > is the 
mean shielding factor. From this formula, one thus sees 
that the product nR can be derived from the measured 
molecular hydrogen fraction /. 

Practically a model is necessary to determine the 
abundance and distribution of H 2 molecules over its ro- 
vibrational levels as a function of depth into the cloud 
and thereby derive the nR product from /. We have used 
an updated version of the stationary model developed by 
Abgrall et al. (1992) and Le Bourlot et al. (1993). The 
model assumes a semi-infinite plane parallel geometry and 
solves the equations of thermal and chemical balance itera- 
tively as a function of depth into the cloud. Transfer in the 
H 2 lines (50 ro- vibrational levels are included in these cal- 
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culations) and dust extinction are taken into account. The 
far UV extinction curves used in the models are specific 
for each star and present a wide variety in terms of 2200 A 
bump and far UV rise. These curves are extrapolations of 
fits to the IUE extinction curves performed by Boulanger 
et al. (1994). The chemical network includes 100 species 
and 775 reactions. To take into account UV penetration 
from two sides, we run models with half the total extinc- 
tion and multiply all integrated column densities by two. 
For all calculations, we have assumed a constant gas den- 
sity through the line of sight. Since we can assume G=l, 
the model has only two free parameters, the gas density n 
and the product nR. The gas kinetic temperature is cal- 
culated in the model at each depth assuming local equi- 
librium between heating and cooling processes. The dom- 
inant heating process is the photo-electric effect on small 
grains. Heating by dissipation of turbulent kinetic energy 
is not included. Model results are shown in Fig. for the 



'. HD108927 
































; n = 20 









nR(s _1 ) 



Fig. 3. Determination of the nR product from the H2 frac- 
tion / for HD 108927. The triangles and square symbols 
represent model calculation for n=20 and 100 cm -3 , re- 
spectively. 

line of sight towards HD 108927, where the computed H 2 
fraction / is plotted versus the product nR. A comparison 
of the model results for n=20 and 100 cm -3 illustrates the 
density independence of the relationship. The horizontal 
lines in the figure represent the measured / with its error 
bar (from Table [l]) . The intersections of these lines with 
the model results define the range for nR. We have used 
the model to quantify the uncertainties on the product 
nR related to the line of sight extinction and the far-UV 
radiation field. We find that the error bar on the Ay mea- 
surement translates into an uncertainty of ±50% on nR. 
An uncertainty of a factor 2 on the far-UV radiation field 
intensity, which is a reasonable assumption, translates into 
an error-bar of the same amplitude on nR. We thus con- 
sider that within the frame work of the model, a factor 
of 2 is the magnitude of the error-bar on nR. Additional 



systematic errors could be coming from the simplifying 
assumptions made in the model. These are of course im- 
possible to quantify. 

The nR values derived for the three stars are listed 
in Table || Their scatter is of a factor of 2, comparable 
to the uncertainty on the individual measurements. Jura 
(1975b) had found a much larger scatter with values from 
510~ 16 to 310~ 14 s -1 . The dispersion of Jura's values 
might be partly due to changes in the G value among the 
lines of sight towards the very luminous stars observed by 
Copernicus. 

3.2. Determination of R from the H 2 density estimate 

To determine the H2 formation rate R one must comple- 
ment the nR values with an estimate of the gas density 
n. 

The gas density n is the second free parameter of the 
model. Density governs the cooling rate and the assump- 
tion of thermal balance makes the gas density and tem- 
perature uniquely related at each depth once the external 
radiation field is given. Any tracer of the gas temperature 
is therefore also a tracer of the density. We thus determine 
the gas density by fitting the model to the column den- 
sity ratio of the two first levels, N(J=1)/N(J=0), known 
to be an indicator of temperature. Fig. |J illustrates the 
best fits, corresponding to the density values indicated in 
each diagram. Note that in the model the temperature is 
not uniform, its value is determined at each depth by the 
thermal balance. The temperature Tb indicated in Fig. ^ 
is that of the shown Boltzmann distribution corresponding 
to the N(J=1)/N(J=0) ratio value. 

The R values derived for the three lines of sight, based 
on these density estimates are listed in Table § This de- 
termination of n assumes that the ortho to para H2 ratio 
is at its equilibrium value and is therefore very sensitive to 
the temperature, which is itself, within the assumption of 
thermal balance, very sensitive to the gas density. Within 
the model framework, the H2 J=l and column density 
ratio constrains density within about 20 %, an uncertainty 
smaller than that on the product nR. The uncertainty on 
R is thus governed by the uncertainty on nR: about a 
factor of 2. 

The values of R for the three lines of sight are close to 
each other and close to the values found by Jura (1975a) 
but with significantly lower uncertainties. Note neverthe- 



Etoile 


HD 102065 


HD 108927 


HD 96675 


nR (s _1 ) 


2.3 10" 15 


0.87 10" 15 


2.0 10" 15 


n(cm~ 3 ) 


50 


28 


50 


^(cmV 1 ) 


4.5 10" 17 


3.1 10" 17 


4.0 10" 17 



Table 3. H2 formation rate R from the product nR and 
the density n estimated from N(H2,J=1)/N(H2,J=0). 



less that the derivation of R from nR is valid within the 
hypothesis that the density is homogeneous in the molecu- 
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Fig. 4. H2 excitation diagrams. The asterisks represent the measured column densities listed in Tablc[T]for J=0 and J=l 
and in Table |^ for the higher J levels within the two assumptions b(H2) C x=b(CH) (higher values) and b(H2) ox =b(CH + ) 
(lower values). The squares represent the model that best fits the J=0 and J=l observations and corresponds to the 
density indicated. The dashed line is the Boltzmann distribution for the temperature Tb corresponding to the observed 
N(J=1)/N(J=0) ratio. 



lar gas. This might not be the case. Indeed, for HD 102065 
for which we have been able to derive three other indepen- 
dent density estimates from the comparison of model cal- 
culations and measured quantities (namely N(Ci), N(CH) 
and the C I fine structure level population from Gry et al. 
(1998)), there is a scatter of a factor of three among the 
four density estimates. This scatter might well reflect a 
density inhomogeneity along the line of sight but could 
also result from other model shortcomings, in particular 
in our understanding of the chemistry in diffuse clouds. 

The nR values directly translate into an H2 forma- 
tion timescale of l/(nR) ~ 2 10 7 yrs. For the three stars 
it is a few times larger than the dynamical timescale 
(~ 10 6 L(pc)/b(km/s)yrs) set by turbulent motions on 
the scales of the absorbing clouds (L = N/n ~ lOpc). 
The model tells us that the photodissociation timescale 
is larger than the H 2 formation timescale especially in 
the shielded layers of the cloud (the dissociation timescale 
goes from about 3 10 7 yrs to about 3 10 8 yrs for Ay from 
0.3 to 1). On the other hand the timescale of ortho to 
para conversion through proton exchange reactions with 
protonated ions (e.g. H + and H^", Gerlich 1990) is much 
shorter (0.2 to 2 10 5 yrs). Consequently, the observed H2 



abundance may not correspond to the equilibrium value 
between formation and destruction and by assuming equi- 
librium in our interpretation we could be overestimating 
R and the nR product. Study of a larger sample of stars 
could reveal scatter in the R and nR values reflecting var- 
ious states of evolution. 

4. Warm H 2 gas 

The excitation diagrams of Fig. [| show that the H2 col- 
umn densities derived from the observations for J > 2 
are all significantly higher than the model values. This 
is true for our two determinations of H2 column densi- 
ties, the highest based on the Doppler parameter of CH 
and the lowest based on that of CH + (see Sect. §. This 
means that H2 excitation by UV pumping and H2 forma- 
tion on grains as computed with the model does not ac- 
count for the J > 2 H 2 column densities. With the model, 
we checked on HD 102065 that to populate these levels by 
UV pumping alone (with no collisional deexcitation) we 
would have to increase the radiation field to a G value of 
17 to reproduce the column densities derived with b(CH) 
and to a G value of 9 to reproduce the column densities 
derived with b(CH + ). Such values are clearly incompati- 
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ble with the IRAS dust data because no enhanced IR dust 
emission is found at the position of the stars and the ra- 
tio between the 100/im cloud brightness and extinction is 
within the range of values observed over the high Galactic 
sky. We have checked also that we cannot reproduce the 
measured J=2 level populations even if all the H2 binding 
energy (4.48 eV) is transformed in H2 internal excitation. 

We thus infer the existence of warm gas along the three 
lines of sight, where the J>2 H2 levels are populated by 
collisional excitation. This statement relies on the validity 
of the H2 column densities determinations. In the absence 
of any direct information, one cannot exclude that the H2 
velocity distribution is broader than that of both CH and 
CH + and consequently that even the lowest set of values 
in Table |^ are larger than the true H2 column densities. 
However, we consider unlikely that the H2 velocity dis- 
tribution could be such that there is no need for warm 
H2 gas because the detection of CH + in the direction of 
both HD 96675 and HD102065 where it has been looked 
for is an independent evidence for the existence of a warm 
H2 component along these lines of sight. If one considers 
only collisional excitation, the H2 excitation temperature 
is lower than the true gas temperature for densities below 
the critical densities of each level. The temperature of the 
warm gas thus needs to be at least equal to the H2 exci- 
tation temperature derived from the J=3 to 5 levels, i.e. 
between 200 and 240 K for the three stars. 

Our conclusion about H2 excitation at J>2 levels dif- 
fers from a common interpretation of UV H2 absorption 
lines (e.g. Jura 1975b) where it is assumed that the high 
J levels are mainly populated by pumping through UV 
photons of the target stars. For our sample of stars we 
are able to rule out this interpretation, because the FUSE 
sensitivity allowed us to observe late B stars which do not 
interact with the absorbing matter, thus lines of sight for 
which the UV radiation field strength is constrained to be 
close to the mean Solar Neighborhood value. 

The existence of warm H2 gas within the diffuse ISM 
has been considered along many lines of sight to account 
for the observed column densities of CH + . The observed 
abundance of CH + is a well known problem of interstel- 
lar chemistry. The only efficient path for CH + formation 
is the highly endothermic (4640 K) reaction between C + 
and H-2. Further, CH + is efficiently destroyed by reaction 
with H2 once it is formed. One thus considers that CH + 
only exists in significant abundance where the molecu- 
lar gas is warm. Away from hot stars, localized volumes 
of warm gas can be created and sustained by dissipation 
of the gas kinetic energy. Formation of CH + have been 
quantitatively investigated in the specific cases where dis- 
sipation occurs within MHD shocks (e.g. Flower & Pincau 
des Forets, 1998) or coherent vortices in MHD turbulence 
(Joulain et al., 1998). In these models the temperature 
of the warm gas and the ratio between CH + and warm 
H2 column densities depends strongly on local physical 
conditions (e.g. shock velocity, gas density, magnetic field 
value). The addition of H2 excitation studies provides a 
mean to constrain these models and can thus help un- 



derstand the physics of kinetic energy dissipation in the 
diffuse ISM. 

For each line of sight we have computed and listed in 
Table |^ a total column density of excited H2 gas (N(H2) e x) 
by summing the column densities in the J=3 to 5 levels. 
The range of values corresponds to the two assumptions 
used for the determination of H2 column densities (see 
Sect. ||). The fact that CH + must form in warm H2 gas 
makes us consider the lowest value of (N(H2) CX , those ob- 
tained for the CH + Doppler parameter, to be the most 
realistic estimate. For HD 108927, in the absence of CH + 
observation only the higher estimate is listed. Note that in 
all cases the excited H2 gas represents a very small fraction 
of the total gas. We have not attempted to separate the 
respective contributions of UV pumping in the cold gas 
and truly warm gas to the column densities of excited H2 . 
The excitation diagrams in Fig. |] show that UV pumping 
is always a minor contribution. 

Based on Copernicus data, Frisch & Jura (1980) sug- 
gested a correlation between the column densities of H2 
in the J=5 level and CH + . They pointed out that CH + 
is particularly abundant in regions with large amounts of 
rotationally excited H2. This result has been confirmed by 
Lambert & Danks (1986) who compared the column den- 
sities N(CH+) and N(H 2 ,J) for the rotational levels J=0- 
5 in a larger sample of sight-lines. They concluded that 
the rotationally excited H2 is a tracer for those hot H2 
molecules that initiate CH + formation. Both lines of sight 
of our study, HD102065 and HD96675, for which the CH+ 
column density is known, fit within the dispersion of the 
Lambert & Danks (1986) correlation diagrams. Note that 
the Lambert & Danks' sample consist of luminous stars 
for which one cannot exclude that the excited H2 gas re- 
sults from the interaction (UV pumping and shock) of the 
star with its surrounding medium. Further sight-lines such 
as those considered here, for which this possibility can be 
excluded, are needed to physically discuss the correlation 
between CH + column densities and warm H2 in relation 
with existing models of CH + formation. 

One Galactic line of sight not crossing any known 
star forming region has been observed with the mid-IR 
spectrometer on board of the Infrared Space Observatory 
(ISO). This observation lead to the detection of three pure 
rotational lines of the H2 ground state (Verstraete et al. 
1999, Falgarone et al. 2002). The H2 populations inferred 
from these detections are compared with our FUSE re- 
sults in Table |[ The N(H 2 ) CX /A V ratio for the long line 
of sight through the Galaxy is a factor 5 and 2 higher that 
the lower estimates listed for HD 96675 and HD 102065. 
We thus find along these two sight-lines a fraction of warm 
H2 gas close to the mean Galactic value derived from the 
ISO observation. The presence of warm H2 thus seems to 
be a general characteristics of the diffuse ISM. 
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HD102065 


HD108927 


HD96675 


ISO 


N(H 2 ) ex 
A v 

N(H 2 )cx 
Av 


1.7- 3.7 10 17 

0.67 
2.5 - 5.4 10 17 


1.0 10 17 

0.68 
1.5 10 17 


0.9 - 1.7 10 17 
1.1 

0.8- 1.6 10 17 


7.7 10 18 

18 
4.3 10 17 



Table 4. Abundance of excited H2 gas. N(H2) CX (in cm 2 ) is the sum of the J=3 to 5 column densities from Table |[ 
The last column (ISO) refers to an infrared observation of a long line of sight through the Galaxy (from Falgarone et 
al. 2002). 



Team by the NASA-CNES-CSA FUSE mission operated by 
the Johns Hopkins University. 
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